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Procede pour evaluer une grandeur physique representative 
d T une interaction entre une onde et un obstacle 

L' invention concerne la modelisation des interactions 
5 entre une onde incidente et un obstacle de cette onde, 
notarament dans le domaine du controle non destructif . 

La modelisation des interactions entre une onde et un 
obstacle recevant cette onde, tel qu'une cible placee dans 
10 la zone sensible d'un capteur, trouve une application 
avantageuse en controle non destructif . 

On connait un procede de modelisation dit "par elements |- 
finis" consistant k appliquer un pavage de l'espace.^ 
15 tridimensional entourant 1 1 obstacle et a evaluer les if 
interactions precitees pour tous les paves de l T espace. ^| 

? - 

Les procedes de calcul par "elements finis" apportent unej, 
resolution a un probleme pose sous forme d' equations:- 

20 dif ferentielles partielles. lis se fondent sur une 
representation de l'espace d' etude par un assemblage 
d f Elements finis, a l'interieur desquels des fonctions 
d' approximations sont determinees en termes de valeurs 
nodales de la grandeur physique recherch6e. Le probleme 

25 physique continu devient done un probldme discret aux 
elements finis ou les grandeurs nodales sont les nouvelles 
inconnues. De tels procedes cherchent done a approcher la 
solution globale, plut6t que les equations aux derivees 
spatiales partielles de depart. 

30 
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La discretisation de l'espace pris en compte assure que ce 
dernier soit entierement recouvert par des elements finis 
(lignes, surfaces ou volumes), cette operation porte le 
nom de "maillage" dans l'espace bidimensionnel (2D) ou de 
"pavage" dans l'espace tridimensionnel (3D). Les elements 
mis en jeu sont soit rectangulaires ou triangulaires en 
2D, soit parallelepipediques ou tetraedriques en 3D. lis 
peuvent etre de tallies differentes, repartis uniformement 
sur la surface ou non. 

En general, la grandeur physique recherchee, telle qu'un 
potentiel electrostatique ou une valeur de pression, est 
connue sur la frontiere du domaine. Cette frontiere peut 
etre fictive. On y impose des conditions aux limites. Le 
potentiel est done inconnu a l'interieur du meme domaine. 
On definit alors un noeud comme etant un sommet d' un 
element. Les inconnues du probleme sont done les valeurs 
du potentiel en chaque noeud de 1' ensemble du domaine. 



A titre d' illustration, la figure 6 de l'art anterieur 
represent e un exemple de surface, constitute par deux 
materiaux Ml et M2, de proprietes electromagnetiques 
differentes, et maillee par des elements triangulaires 
comportant chacun trois nceuds Ai, Bi et Ci. L' ensemble du 
25 domaine est delimite par une frontiere F. 



Une fois le maillage defini, plusieurs approches existent 
pour transformer la formulation physique du probleme en 
une modelisation discrete aux elements finis. Si le 
30 probleme est formule par des equations dif f erentielles et 
consiste a minimiser une f onctionnelle, on applique 
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gen§ralement une methode variationnelle- Cette 
transformation aboutit a une formulation matricielle dont 
. la resolution donne les solutions nodales, les solutions 
aux points non nodaux etant obtenues par interpolation 
5 lineaire, 

Neanmoins, de tels calculs, en trois dimensions, 
requi^rent des ressources inf ormatiques iraportantes et 
generent des temps de calcul tres longs , en depit de 
10 1 1 accroissement des performances de logiciels permettant 
d' implement er ces calculs. 

Certes, les problemes 2D, souvent simplifies par des^ 
conditions de symetrie avantageuses pour ne modelisez^ 

15 qu'une partie de la geom<§trie, se resolvent rapiderrient 

Mais il n'en est rien pour les problemes 3D, qui sont les ; j| 
plus frequents. La figure 6 montre combien la finesse du^ 
maillage, c'est a dire le rapport entre la taille d'ur£. 
element et celle du plus petit detail du domaine, fait 

20 • croitre le nombre de noeuds. 

En consequence, la quantity d' equations et . d'inconnues 
augmente proportionnellement , et, de la, le temps de 
calcul necessaire a la resolution du probleme. II est 
25 important de pr6ciser que la generation du maillage, & 
savoir la discretisation de l'espace de travail, et la 
generation de la liste des noeuds consomme un temps de 
calcul superieur & celui necessaire a la resolution du 
probleme . 

La presente invention vient am^liorer la situation. 
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Elle propose a cet effet un precede pour evaluer une 
grandeur physique associee a une interaction entre une 
onde et un obstacle, dans une region de l'espace 
tridimensionnel, dans lequel : 

a) on determine, par maillage, une pluralite 
d'echantillons de surface dont une partie au moins 
represent e la surface d'un obstacle recevant une onde 
principale et emettant, en reponse, une onde secondaire, 
et l'on attribue a chaque echantillon de surface au moins 
une source emettant une onde elementaire representant une 
contribution a ladite onde secondaire, 

b) on forme un systeme matriciel comport ant : 

- une matrice d' interaction, inversible, appliquee a 
une region donnee de l'espace et comportant un nombre 
de colonnes correspondant a un nombre total de 
sources, 

- une premiere matrice colonne dont chaque coefficient 
est associe a une source et caracterise l'onde 
elementaire qu'elle emet, 

- et une seconde matrice colonne, obtenue par une 
multiplication de la premiere matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d'une grandeur physique representative de 
l'onde emise par 1' ensemble des sources en ladite 
region donnee, 

c) pour estimer les coefficients de la premiere matrice 
colonne, on affecte des valeurs de grandeur physique 
choisies a des points predetermines, associes chacun a un 
echantillon de surface, lesdites valeurs choisies etant 
rangees dans la seconde matrice colonne, et l'on multiplie 
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cette seconde matrice colonne par 1' inverse de la matrice 
d 1 interaction appliquee auxdits points predetermines, 
d) pour evaluer ladite grandeur physique repr^sentant 
l'onde emise par I 1 ensemble des sources dans une region 
5 donn6e de I'espace tridimensionnel, on applique la matrice 
d' interaction a ladite region donnee et on multiplie cette 
matrice d' interaction par la premiere matrice colonne 
comportant les coefficients estimes a 1'etape c) • 

10 Ainsi, selon I'un des avantages que procure la presente 
.invention, 1'etape de maillage a) ne concerne qu'une ou 
plusieurs surfaces, tandis que le procede de modelisation^ 
de type "par Elements finis" necessite un pavage de tout^ 
I'espace avoisinant 1" obstacle, ce qui permet de reduire,^ 

15 dans la mise en ceuvre du procede selon l 1 invention, les S. 
ressources memoires et les temps de calcul necessaires. ^ 

y 

Le procede selon l 1 invention s' applique aussi . bien & une.£ 
onde principale emise par une source lointaine qu'a une 
20 onde principale emise en champ proche. 

Avantageusement, pour evaluer une grandeur physique 
representative d'une interaction entre un element 
rayonnant une onde principale en champ proche et un 

25 obstacle recevant cette onde principale, 

- £ 1'etape a), on determine en outre, par maillage, une 
pluralite d 1 6chantillons de surface representant ensemble 
une surface active de l'616ment rayonnant I'onde 
principale, et l'on attribue & chaque echantillon de la 

30 surface active au moins une source emettant une onde 



6 

elementaire representant une contribution a ladite onde 
principale, 

- on applique les etapes b) , c) et d) aux echantillons de 
la surface active, et 

5 - on evalue ladite grandeur physique representant 
l f interaction entre 1' element rayonnant et 1 T obstacle dans 
une region donnee de l'espace tridimensionnel, en prenant 
en cornpte la contribution, en ladite region donnee, de 
l'onde principale emise par 1' ensemble des sources de la 

10 surface active et la contribution de l'onde secondaire 
emise par l 1 ensemble des sources de la surface de 
1' obstacle. 

On entend par les termes "element rayonnant" aussi bien un 
15 emetteur de l'onde principale, tel qu T un g^nerateur 
d'ondes, qu'un r6cepteur de l'onde principale, tel qu'un 
capteur de cette onde. 

Dans un premier mode de realisation, la grandeur physique 
20 a evaluer est une grandeur scalaire et, a l'etape a), on 
attribue une source unique a chaque echantillon de 
surface. 

Dans un second mode de realisation, la grandeur physique a 
25 6valuer est une grandeur vectorielle exprim^e par ses 
trois coordonndes dans l'espace tridimensionnel, et on 
attribue, a l'etape a), trois sources a chaque echantillon 
de surface. 

30 Dans une realisation avantageuse, pour estimer, a l'etape 
d) , la contribution de l'onde secondaire dans la region 



donnee de . l'espace, les valeurs de grandeur physique 
choisies a l'§tape c) sont fonction d'un coefficient 
predetermine de reflexion et/ou de transmission de I'onde 
principale par chaque echantillon de surface de 
obstacle. 

Ainsi, on comprendra que I'onde secondaire peut aussi bien 
correspondre a une reflexion de I'onde principale, qu'a 
une transmission de I'onde principale, ou encore a une 
diffraction de I'onde principale. Dans cette realisation 
avantageuse, I'etape c) correspond finalement & une 
determination des conditions aux limites a la surface de 
I 1 obstacle, en tant qu' interface entre deux milieux 
distincts notamment dans une heterostructure . 

Par ailleurs, pour un controle non destructif d T une cible 
formant un obstacle d'une onde principale, on attribue un 
coefficient de reflexion ou de transmission choisi a tous 
les points predetermines de la surface de la cible, et on 
compare une . simulation obtenue par la mise en ceuvre du 
procede au sens de 1' invention avec une mesure 
exp6rimentale. Ainsi, les points de la surface de la cible 
qui, dans la mesure experimentale, ne v^rifient pas la 
simulation correspondent a des inhomog6n<§ites ou a des 
impuretes sur la surface de la cible. 

Dans une autre approche, on connait les propri£t6s 
globales de 1' obstacle, notamment en transmission et/ou en 
reflexion. Par la mise en ceuvre du proced6 de 1' invention, 
on optimise alors la position dans I'espace d'un capteur 
ou m§me la forme de ce capteur, pour une application a un 
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controle non destructif, ce capteur etant destine a 
analyser une cible formant un obstacle de l'onde 
principale. 

A cet effet, dans une realisation avantageuse, on compare 
une pluralite de valeurs de la grandeur physique estimee a 
1'etape d) du precede au sens de 1' invention, obtenues 
pour une pluralite de regions de l'espace, pour 
select ionner une region candidate pour la disposition d'un 
Element rayonnant destine a interagir avec 1' obstacle. 

Comme indique ci-avant, on entend par les termes "element 
rayonnant" aussi bien un capteur qu'un generateur de 
l'onde. On comprendra ainsi que 1 ' optimisation de la 
position de 1' element rayonnant peut s'appliquer aussi a 
1' optimisation de la disposition ou de la forme d'un 
generateur d'ondes. Par exemple, la presente invention 
trouve encore une application avantageuse a la disposition 
de haut-parleurs dans un volume ferme, delimite par des 
obstacles, oomme par exemple 1'habitacle d'un vehicule 
automobile- 

D'autres caracteristiques et avantages de 1' invention 
apparaitront a l'examen de la description detaillee ci- 
apres, et des dessins annexes sur lesquels : 
- la figure 1A represente schematiquement les surfaces 
respectives d'un el6ment rayonnant ER emettant une onde 
et d'un obstacle OBS recevant cette onde, maillees en 
vue d'evaluer une grandeur scalaire representative de 
l'onde en un point M de l'espace tridimensionnel ; 
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- la figure IB represente en detail un echantillon de 
surface dSi correspondant a une maille de la figure 1A, 
ainsi qu'une source Si associee a 1 T echantillon de 
surface dSi ; 

5 la figure 2A represente schematiquement les surfaces 

respectives d ? un element rayonnant ER ernettant une onde 
et d'un obstacle OBS recevant cette onde, maillees en 
vue d'evaluer une grandeur vectorielle representative 
de l'onde en un point M de l T espace tridimensionnel ; 

10 - la figure 2B represente en detail . un echantillon de 
surface dSi correspondant a une maille de la figure 2A, 
ainsi que trois sources associees SAi, SBi et SC± ; 

- la figure 2C represente, en vue de face, une surface*,, 
maiilee dont chaque echantillon de surface comportef 

15 trois: sources SA±,. SBi et SCi, pour 1' estimation d'une| 

grandeur vectorielle ; ~ f£ 

- la figure 3A represente, a titre illustratif, lesp 
armatures d'un condensateur , de potent iels elect riques? 
respectifs VI et V2, pour 1 1 estimation d'un potentiel 

20 electrique au point M de l'espace tridimensionnel, a 

chaque echantillon de surface dSi de la figure 3A etant 
associee une unique source Si ; 
*- la figure 3B represente, a titre illustratif, les 
armatures d'un condensateur, de potentiels eiectriques 

25 respectifs VI et V2, " pour l'estimation d'un champ 

electrique E(M) , au point M de l'espace 
tridimensionnel, a chaque echantillon de surface dSi de 
la figure 3B etant associee trois sources SA if SB ± et 
SCi ; 
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- la figure 4 A represente, comme les figures 1A et 2A, 
une interaction entre un element rayonnant ER et un 
obstacle OBS, pour evaluer une grandeur physique 
(scalaire ou vectorielle) en un point M dans une 
5 portion de l'espace delimitee par la surface de 

1' element rayonnant et la surface de 1' obstacle, ce 
point de l'espace M recevant a la fois I'onde emise par 
l 1 element rayonnant et I'onde reflechie par 
1 ' obstacle ; 

10 - la figure * 4B, complementaire de la figure 4 A, 
represente une transmission par 1' obstacle OBS de 
l f onde emise par 1' Element rayonnant ER, en un point M 
d'un demi-espace delimite par le plan que forme la 
surface de l f obstacle OBS ; 

15 - la figure 5A represente sch§matiquement un obstacle 
OBS, de dimensions finies, avec des sources associees 
aux echantillons de surface agencies pour estimer une 
grandeur representative d f une reflexion de I'onde sur 
1' obstacle ; 

20 - la figure 5B, en complement de la figure 5A, represente 
schematiqueraent un obstacle OBS, de dimensions finies, 
avec les sources associees aux Echantillons de surface 
disposees pour estimer une grandeur representative de 
la transmission de I'onde par 1' obstacle ; 

25 - la figure 5C represente une simulation d'une onde 
ultrasonore emise par un element rayonnant ER et se 
propageant vers un obstacle OBS ; et 
- la figure 6 represente un maillage de milieux 
tridimensionnels, pour 1 T application d'un procede de 

30 calcul par "elements finis", au sens de l'6tat de la 

technique. 
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On se refere tout d'abord a la figure 1A, sur laquelle la 
surface d'un obstacle OBS, recevant une onde, est maillee 
selon une pluralite d f echantillons de surface dSi a dS 4 , 
5 conf ormement a l'6tape a) precit6e. 

En se referant a la figure IB, a chaque echantillon de 
surface dSi est associe un hemisphere HEMi, tangent a. 
1 1 echantillon de surface dS± en un point de contact 

10 Pref 6rentiellement , ce point de contact Pi correspond au 
soinmet de 1 1 hemisphere HEMi . Pour l 1 estimation d'une 
grandeur physique scalaire au point M (tels qu T un ^ 
potentiel electrostatique, une pression acoustique ou ^ 
autre) , une source unique Si est associee & 1 1 echantillon i 

15 ' de surface dSi. Comme on le verra plus loin, dans le cas 
de 1' estimation d'une grandeur vectorielle dans un point 
de l'espace M, on affectera plutot trois sources a chaque ^ 
echantillon de surface dS±. 

20 * Pref erentiellement , 1 T hemisphere HEMi est construit comme 
decrit ci-apr&s. Pendant l'etape de maillage a) precitee, 
on lvalue, d'une part, la surface de l 1 obstacle OBS et, 
d" autre part, on choisit un nombre d' echantillons de 
surface dSi selon la precision souhaitee de 1' estimation 

25 de la grandeur physique au point M. Ainsi, la surface d'un 
Echantillon dSi est donnee par S 0 /N, oft S Q correspond £ la 
surface totale de 1' obstacle et N correspond au nombre 
choisi d' echantillons de surface dS;*.. 
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L ' hemisphere HEMi est de meme surface que 1 ' echantillon 
dSi. Ainsi, le rayon Ri de 1' hemisphere se deduit de 
l f expression : 



Chaque maille que represente un echantillon de surface d Si 
presente, dans 1- example decrit, une forme de 
parallelogramme, de centre Pi correspondant au point 
d' intersection des diagonales de ce parallelogramme. 
L' hemisphere HEMi est tangent a 1 • echantillon de surface 
dSi en ce point P*. Bien entendu, les mailles peuvent etre 
de forme differente, triangulaire ou autre. On indique de 
facon generale que le point Pi correspond au barycentre de 
la maille. 

On definit ainsi la position de la source Si (situee au 
centre de 1' hemisphere HEMi) . La distance separant la 
source S, du point de contact Pi correspond au rayon Ri de 
1< hemisphere HEMi et la droite qui passe par les points Pi 
et Si est orthogonale a la maille dSi- 

Dans 1- example represente sur la figure 1A, on maille en 
outre la surface d'un element rayonnant BR, correspondant 
par exemple a un generateur d'ondes. A chaque maille de la 
surface de 1' element rayonnant est associe un echantillon 
de surface dS'i, comme on le verra plus loin. 

Le systeme matriciel que l'on met en forme a 1'etape b> 
precitee correspond a : 
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'vcmo > 




(y x \ 


V(M 2 ) 


= F x 


v 2 


,V(M N ), 







[1] 



oh : 

les coefficients Vj (avec j— 1,2,...,N) de la premiere 
matrice colonne correspondent a des valeurs associees 
chacune a une source, telle qu'une charge electrique 
(dans le cas de 1' estimation d'un potentiel 
electrique) , ou encore k un flux magnetique (dans le 
cas de l'estimation d T un potentiel magnetique), ou : 
encore a une vitesse sonore (dans le cas de 
l'estimation d'une pression acoustique liee & la 
propagation d'une onde sonore) ; 
- les coefficients V(M ± ) (avec i=l,2,...,N) de la seconde 
matrice colonne correspondent chacun a une valeur de la. 
grandeur physique (un potentiel electrique ou 
• magnetique ou une pression) £ estimer en un point Mi de* 
I'espace ; 

la matrice d' interaction F comporte des coefficients 
Ci,j dont 1' expression gen§rale est donnee par : 

Gj=f(M,Sj) < [2] 

On comprendra ainsi que les coefficients de la matrice F 
sont des coefficients d f interaction qui dependent de la 
distance s^parant chaque point de I'espace M± d'une source 
Sj associee a une maille dSj. 
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Dans le cas de la propagation d'une onde electrique, lea 
coefficients c M , v 3 et V( Mi >, respectivement de la 
matrice d' interaction de la premiere et de la seconde 
matrice colonne, sont donnes par : 



1 

C- • = = 

1,J 2xe 0 M i S j 

[3] 

vj =qj 



V(Mj) = Ui 



ou : 



- so est une constante dielectrique, 

- MjSj est une distance mesuree en valeur algebrique, 

- correspond a une charge electrique caracterisant une 
source Sj, et 

- Ui correspond a un potentiel electrique au point M*. 

Dans le cas de la propagation d'une onde magnetique, 
1' expression de ces coefficients est la suivante : 



Ci »- 



ou 



1,J 2^o M i S j 

[4] 

Vj =<p } 
V(M L ) = 0j 

m correspond a la permeabilite magnetique du milieu ou 
se situe le point Mi, 
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- q>j correspond au flux magnetique associe a la source 
Sj ; 

- 0i correspond au potentiel magnetique au point Mi. 

Dans le cadre de la propagation d'une onde ultrasonore, 
ces coefficients sont donnes par : 

- " 15] 

v j =v j 

V(Mi) = Pi 



10 dans laquelle : 

X — X . * 

- a> est la pulsation de l T onde sonore ; f. 

- p est la density du milieu dans lequel se situe le 
point Mi ; 

15 - le vecteur Vj correspond £ la vitesse sonore issue de 
la source Sj ; 

- k correspond au vecteur d'onde de l'onde sonore ; et 

- pi correspond a la pression acoustique gen^ree par la 
propagation de l'onde ultrasonore au point Mi- 



Dans 1 T expression des coefficients ci,j, le terme dSj 
correspond a la surface de 1 1 6chantillon associe a la 
source S j • De preference, le maillage d'une surface au 
sens de l'etape a) du procede selon 1" invention, est 
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choisi de sorte que chaque maille comporte une merne 
surface dS=dSi=dS 2 =...=dSj . 

On remarque en particulier dans 1' expression des 
coefficients Ci fj qu'ils dependent du produit scalaire 
entre le vecteur d'onde et le vecteur M^j . Ainsi, pour 
des ondes ultrasonores, on tient compte d'un dephasage 
entre les chemins qui lient chaque source Sj a un point de 
l'espace tridimensionnel M, ce dephasage etant du a une 
difference de marche entre les rayons partant de chaque 
source et arrivant au point M (comrae le raontre la figure 
4A) . En particulier, on comprendra que 1' angle d 1 incidence 
d'un tel rayon est pris en compte dans 1' expression des 
coefficients de la matrice d' interaction F. 

Bien entendu, dans le cadre de la propagation d'une onde 
electromagnetique de frequence elevee, done de longueur 
d'onde courte, qui differe du cadre electrostatique ou 
magnetostatique ci-avant, on tient compte du terme de 
propagation expjjf£?) dans 1' expression de la matrice 
d' interaction, par rapport a la geometrie du probleme a 
resoudre, comme dans le cadre de la propagation d'une onde 
ultrasonore ci-avant (relation [5] ) - 

Ainsi, le syst^me matriciel de 1' equation [1] permet 
d'estimer, a partir d'une matrice d' interaction F et d'un 
vecteur comportant les valeurs v 3 associees aux sources 
Sj, les coefficients d'un vecteur (matrice colonne) 
comportant les valeurs de grandeur physique V(Mi) aux 
points de l'espace Mi- 
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Pour determiner les valeurs des sources Vj, on applique le 
systeme matriciel suivant : 







f Vl l 


V(P 2 ) 


= F x 




,V(P N ), 







[6] 



ou : 

les coefficients de la matrice d' interaction F 
s'expriment par avec Gj—f(BSj) 

les indices i et j correspondent respectivement £ la 
i toe ligne et la j 6me colonne de la matrice d* interaction 
F. Cette matrice d 1 interaction , comporte, pour la 
determination des valeurs associees aux sources Vj, N 
lignes et N colonnes, en rappelant que N est le norabre 
total de mailles a la surface de 1' obstacle ; 
- les points Pi correspondent au sommet des hemispheres 
HEMi de la figure IB. 

La raise en ceuvre de l f <§tape c) du procede au sens de la 
presente invention correspond a calculer une condition aux 
limit es pour les points P ir de proprietes connues, comme 
on le verra plus loin. 

Le systeme matriciel de 1' equation [6] devient alors : 
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'V(P,) > 


v 2 


= F -1 x 


V(P 2 ) 

• 






lv(P N ), 



[7] 



ou : 

- F" 1 correspond a 1' inverse de la matrice d' interaction 
F ; et 

- les valeurs V(Pi) sent predetermines, en fonction des 
conditions aux limites precitees. 

On determine ainsi les valeurs de source Vj . 

A partir de 1' estimation de ces valeurs de source Vj , on 
peut calculer la grandeur physique scalaire en un point M 
quelconque de l'espace tridimensional, a partir de la 
relation : 



V(M)=f]f(MSj)vj 



[8] 



Pour obtenir cette expression de la grandeur scalaire 
V(M), la matrice d' interaction F peut ne comporter qu'une 
ligne de coefficients c-j, avec : 

Cj = f(MSj) , 
mais comporte toujours N colonnes. 

En se referant a nouveau a la figure 1A, on comprendra que 
la surface de 1' obstacle OBS recevant l'onde qu'emet 
!• element rayonnant ER agit, elle-meme, comme une surface 
active re-emettant une onde secondaire (par exemple par 



19 



reflexion) . Chaque source Si represente une contribution a 
1' Emission de cette onde secondaire. 

En outre, pour tenir compte a la fois de la presence de 
l'onde principale et.de la presence de l'onde secondaire 
au point M, on estime la contribution de l'onde principale 
et la contribution de 1 * onde secondaire au point M par le 
systeme matriciel : 



10 



15 



20 



25 



V(M) = Fx 



hi 






v 2 


+ Fx 











[9] 



ou : 

- F» est la matrice d' interaction entre la surface de*, 
1' Element rayonnant et le point M ; 

- v'j (avec j=l,2,3,...,N' ) est la valeur des sources 1 
attributes a chaque tchantillon de surface dS'j de 
l'§16ment rayonnant, N' etant le nombre total de 
mailles choisi pour la surface active de 1' element 
rayonnant ER. 

Les coefficients de la matrice F' sont encore fonction de 
la distance MS'j, ou S*j sont les sources affectees a 
chaque echantillon de surface dS'j de 1' element rayonnant. 

Selon une caracteristique . avantageuse, les valeurs des 
sources de 1* obstacle vj sont determinees en fonction des 
valeurs des sources de 1' element rayonnant v'j, lesquelles 
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sont elles-memes calculees conime on le verra plus loin en 
reference aux figures 4A, 4B, 5A et 5B. 

On se refere maintenant a la figure 2A, dans laquelle on 
5 affecte trois sources & chaque echantillon cie surface dSi, 

en vue d'estimer une grandeur physique vectorielle V(M) , 
en un point M de l'espace tridimensionnel • 

En effet, on comprendra que pour estimer la grandeur 
10 vectorielle, par ses trois coordonn<Ses dans l'espace x, y 
et z, le nombre d T equations a resoudre par rapport au 
systeme matriciel precedent doit etre multiplie par trois, 
Ainsi, la matrice F" 1 de la relation [7] doit comporter 
trois fois plus de lignes que pr6cedemment . La matrice 
15 d 1 interaction F doit, elle-m§me, comporter trois fois plus 
de colonnes que pr^cedemment et 7 a cet effet, on prtvoit 
avantageuseraent trois sources par maille lorsqu'il s T agit 
de determiner les coordonnees dans l f espace 

tridimensionnel d'un vecteur V(M) • 

20 

En se r6f6rant a la figure 2B, les trois sources SA±, SBj., 
SC lf attributes k un echantillon de surface dSi sont de 
positions respectives determinees comme on l T indique ci- 
apres. Telles que representees sur la figure 2B, les trois 
25 sources SA if SB if SCi sont coplanaires et le plan 
comportant ces trois sources comporte en outre la base de 
1' hemisphere HEMj., L 1 hemisphere HEMi est construit comme 
indique ci-avant (de m£me surface que la surface de la 
maille), avec toutefois le centre de l f hemisphere qui 
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correspond ici au barycentre des trois sources SA±, SBi et 

On entend par "centre de 1 'hemisphere" le centre du disque 
5 qui constitue la base de 1 1 hemisphere . 

Les trois sources qui sont attributes a 1 1 echantillon de 
surface dSi sont disposees aux sorranets d'un triangle 
equilateral dont le barycentre Gi correspond au centre de 
10 1 'hemisphere. De preference, chaque source SAi, SBi et SCi 
est disposee sur le milieu d'un rayon Ri de 1 "hemisphere. 
Ainsi, les droites qui relient le barycentre Gi & chaque 
source sont 6cartees angulairement de 120°. ."f 

15 En se reftrant a la figure 2C, 1 f orientation angulaire des* 
triangles que forment les triplets de sources est choisie y 
aleatoirement, d'un Echantillon de surface a l" autre..'/ 
Avantageusement, on evite ainsi des artefacts de sur- : 
periodicite, dans 1' estimation de la grandeur vectorielle 

20 au point M, qui pourraient resulter du choix d'une meme 
orientation angulaire de ces triangles. 

En reference avec les different s types d'ondes indiques 
precEderament, la grandeur vectorielle V(M) a estimer peut 
25 §tre : 

un champ electrique, dans le cadre de la propagation 
d'une onde Electrique ; 
- un champ magnet ique, dans le cadre de la propagation 
d'une onde magnet ique ; et 
30 - une vitesse du son au point M, dans le cadre de la 
propagation d'ondes ultrasonores . 
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Pour determiner lea valeurs associees a chaque source SAi, 
SBi, SC±, on met en forme le systeme matriciel selon la 
relation suivante : 







vA, > 


v x (P 2 ) 




vA 2 


V X (P N ) 




vA N 


Vy(Pl) 






V y (P 2 ) 




vB 2 


Vy(P N ) 




vB N 


V z (Pi) 




vC, 


V Z (P 2 ) 




vC 2 


,V Z (P N ). 


> 





[11] 



On remarque, en particulier, que la matrice d' interaction 

est de dimensions 3Nx3N, ou N est le nombre total 

d'echantillons de surface. La matrice d' interaction 
s'exprime ici par la relation : 



N 



N{c A (i,j) C|(i,j) Cfc(i,j) 
N^(i,j) Cj&j) C£(i,j) 
N{c A (i,j) Cl(i,j) C£(iJ) 



[12] 
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Les coefficients de cette matrice s ' expriment par : 

CZ{iJ) = fuWPi,Saj)} [13] 

5 avec or = A, B, C 

i = 1, 2, N 

j - 1, 2, ... r N 

u = x, y, z. 

10 et sont encore fonction d'une distance separant le point 
de contact Pi de l T une des sources SOj (cr « A, B ou C) 
d'un triplet associe a un echantillon de surface dSj. 

Par inversion de la matrice d' interaction Fy , on 
15 determine ainsi les valeurs vOj associ6es a chaque source 
SOj, en appliquant des conditions aux limites sur les 
valeurs du vecteur V aux points P±. Ces conditions aux. 
limit es imposent une valeur du vecteur V , selon ses trois 
coordonn6es V x (Pi), V y (Pi) et V 2 (Pi). 

20 

Une fois que ces valeurs de sources vc^ sont ainsi 
determines, on calcule facilement I'expression du vecteur 
V en un point quelconque M de l'espace, par la relation : 

25 V(M) = V x (M)x + V y (M)y + V z (M)z [14] 

W= T,fx[{d{M 9 Saj)}.vaj 
cr=A,B,C 
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j=\...,N 
cr—A^ByC 

j=l,...,N 
<r=A,B,C 



5 x, y et z correspondent a des vecteurs unitaires portes 
par les axes x, y et z de l'espace tridimensionnel . 

Ainsi, la matrice d' interaction ?y , lorsqu'elle est 
appliquee a un point quelconque M de l'espace, ne comporte 
10 finalement que trois lignes associees chacune a une 
coordonnee de l'espace x, y ou z. 

Pour differents types d'ondes, les valeurs des sources V05 
sont, comme precedemment, une charge electrique pour une 
15 onde electrique, un flux magnet ique pour une onde 
magnetique, une vitesse du son pour une onde ultrasonore. 

Plus precisement, les coefficients de la matrice 
d' interaction se determinent a partir des relations 

20 [3], [4] et [5] precedentes en precisant toutefois que : 

V(M) = -^d[V(M)] [15] 



V(M) etant la grandeur scalaire calculee precedemment par 
25 1' equation [8] . 
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Ainsi, pour 1« estimation d'une grandeur vectorielle V au 
point M et pour les types d'onde cites ci-avant a titre 
d'exemple (electrique, magnetique et ultrasonore) , les 
coefficients de la matrice d' interaction Fy sont 
inversement proportionnels . au carre d'une distance 
separant chaque source du point M, alors que pour 
1' estimation d'une grandeur scalaire V en un point M de 
l'espace, les coefficients de la matrice d* interaction F 
sont simplement inversement proportionnels a cette 
distance. Chaque distance implique l'une des sources d'un 
triplet d'un echantillon de surface et un point M de 
l'espace. La matrice d' interaction F y comporte alors 3N^ 
colonnes lorsqu'il s'agit de prendre trois sources par| 
echantillon de surface, alors que la matrice d' interaction^ 
F pour 1' estimation de la grandeur scalaire ne comportait^ 
que N colonnes puisque seule une source par echantillon de^ 
surface etait necessaire. •„.-;. 

I; 

Plus generalement, on attribue une source par echantillon 
lorsque l'on connait des conditions aux limites pour une 
grandeur scalaire et trois sources par echantillon lorsque 
l'on connait plutdt des conditions aux limites pour une 
grandeur vectorielle. 

On se refere maintenant a la figure 3A pour decrire, a 
titre illustratif, une application du procede selon 
1" invention a 1 ' estimation d'un potentiel electrique en un 
"point M de l'espace tridimensionnel, situe entre deux 
armatures d'un condensateur . Les armatures de ce 
condensateur sont portees a des potentiels respectifs VI 
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et V2. La mise en oeuvre de l f 6tape a) consiste d'abord a 
mailler les surfaces respectives des deux armatures* Dans 
l'exemple repr§sente sur la figure 3A, on n'a represents 
que deux mailles pour chaque armature, simplement a titre 
5 illustratif • 

L' application de l'etape b) consiste a mettre en forme le 
systeme matriciel impliquant la matrice d' interaction F et 
le vecteur colonne comportant les valeurs des sources Si & 
10 S 4 . La multiplication de ces deux matrices permet 
d'obtenir un vecteur colonne comportant les valeurs du 
potentiel en un ou plusieurs points M de l'espace. 

La mise en ceuvre de l f 6tape c) du procede selon 
15 1' invention consiste & appliquer le systeme matriciel aux 
points de contact des hemispheres Pi a P 4/ de chaque 
echantillon de surface dSi a dS 4 . H en r^sulte la 
relation suivante : 







f 1 


1 


1 


1 ^1 








PlSi 

1 


P1S2 

1 


P1S3 

1 


P1S4 

1 


hi 


v(p 2 ) 


1 


P 2 Si 


P2S2 


P 2 S 3 


P2S4 


v 2 


v(p 3 ) 
Iv(p 4 )J 


2ft£Q 


1 

P3S1 

1 


1 

P3S21 

1 


1 

P3S3 

1 


1 

P3S4 

1 


v 3 

< v 4; 






IP4S1 


P4S2 


P4S3 


P4S4, 





[16] 



avec V(P!) = V(P 2 ) = Vi 

V(P 3 ) = V(P 4 ) - V 2 
Vi = qi, v 2 = q 2 , v 3 = q 3 / v 4 = q 4 
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Ici, la condition aux limites impose que la valeur du 
potentiel aux points de contact Pi et P 2 correspond au 
potentiel VI de la premiere armature. De meme, le 
potentiel electrique aux points de contact P 3 et P 4 
correspond au potentiel electrique de la seconde armature 
V2. Par inversion de la matrice d' interaction appliquee au 
point de contact Pi, on determine les valeurs des sources 
vi qui correspondent, comme exprime dans la relation [16], 
a des charges 61ectriques qi. 

Les coefficients de ia matrice d' interaction 2yrg * p . g , sont 

parfaitement connus, puisque les positions des sources S j; ; 
et. les positions des points de contact Pi sont determinees^ 
■au prealable, comme represente sur la figure IB. . ^ 

15 - • ...... .i'j 

L ' expression du potentiel electrique V(M) au point M entre^ 
les deux armatures est finalement donnee par_. 
1 ' expression : 
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V(M) = _!_f-2L_ + JBj2_ + _i3_ + -^4_ > | [17] 
K J 2^o l MS l MS 2 MS 3 MS 4j 



On se refere maintenant a la figure 3B sur laquelle on a 
represente les memes armatures que sur la figure 3A, avec 
sensiblement un meme maillage, mais dans le but, ici, 
25 d'estimer une grandeur vectorielle correspondant au champ 
electrique E(M) , au point M de l'.espace tridimensionnel . 
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On peut appliquer les relations [11] a [15] pour estimer 
la valeur du champ electrique au point M, avec, dans la 
relation [13] : 



2ns Q d 2 (P h S<Tj) 



[18] 
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avec 



a = A, B, C 

u = x, y, z 

i = 1, 2, 3, 4 

j = 1, 2, 3, 4, 



Toutefois, les valeurs du champ Electrique au point de 
contact Pi restent a determiner dans la relation [11] . 



15 On introduit alors une loi generale predeterminee du 
comportement du champ (en reflexion, en transmission, ou 
autre) au niveau de la surface de 1' obstacle (des 
armatures dans l'exemple du condensateur precite) , pour 
connaitre les valeurs des sources Vo-j. 

20 

Par exemple, si l'onde electrique est totalement refl^chie 
par la surface d'un obstacle (par exemple l'une des deux 
armatures), le champ electrique en un point de contact P± 
est normal a la surface dSi et ses composantes E x et E y 
25 sont nulles. A titre illustratif, si la surface de 
l 1 armature n'etait representee que par un seul echantillon 
de surface a trois sources, les valeurs de ses sources vA, 
vB et vC seraient toutes egales entre elles a une meme 
valeur +q. 
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Au contraire, si le coefficient de reflexion est 
pratiquement nul a la surface dSi, la composante du champ 
elect rique E 2 au point Pi est nulle, ce qui correspond 
bien au cas ou le champ est sensiblement tangent a la 
surface dSi- Ainsi, a titre illustratif, si la surface de 
1' armature n'etait representee que par un seul echantillon 
de surface a trois sources, les valeurs de ses sources vA, 
vB et vC seraient respectivement, par exemple, +q, +q et 
-2q. Par exemple, dans le cadre de la propagation d'une 
onde magnetique, si la surface d'un capteur a courants de 
Foucault (avec une composante normale du champ magnetique. 
nulle) etait representee par un unique echantillon de 
surface, les flux magnetiques des trois sources associees^ 
& cet echantillon de surface seraient -Kp,+9et -2q>. £ 

On comprend ainsi qu'avec les trois sources parj, 
echantillon dSi, on peut definir, par exemple en fonction, 
de la ponderation de chaque source, une orientation 
quelconque du champ a la surface de 1' obstacle. 

Bien entendu, cette demarche suppose que le coefficient de 
reflexion R d'un obstacle soit connu auparavant. En 
particulier, il peut etre avantageux de comparer une 
simulation et une mesure experimentale pour detecter, a la 
surface d'un obstacle, des inhomogeneites ou des impure tes 
qui correspondent a des points de la surface de cet 
obstacle qui ne verifient pas les valeurs du coefficient 
de reflexion R imposees a chaque point predetermine Pi de 
1 ' obstacle . 
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On peut affecter ainsi une valeur predeterminee du 
coefficient de reflexion a chaque point Pi de la surface 
de l 1 obstacle. A cet effet, on introduit une matrice R qui 
est representative du coefficient de reflexion en chaque 
5 point Pi. Pour une interaction entre un element rayonnant 
et un obstacle , on peut ainsi exprimer le systeme 
matriciel de la relation [9] autrement, c'est-a-dire en 
donnant une expression unique de toutes les sources du 
systeme (a la fois de l 1 obstacle et de 1" element 
* 10 rayonnant), comme indique ci-aprds. 

Dans ce qui suit, on indique que : 
• - F(P) est la matrice d T interaction de l 1 obstacle OBS 
appliqu6e aux points Pi de la surface de 1' obstacle 
15 OBS ; 

- F(P') est la matrice d' interaction de I 1 obstacle OBS 
appliquee aux points P'i de la surface de 1' Element 
rayonnant ER; 

- F 1 (P) correspond a la matrice d f interaction de 
20 l'616ment rayonnant ER appliquee aux points P ± de la 

surface de 1* obstacle OBS ; 

- F' (P') correspond a la matrice d T interaction de 
1' element rayonnant ER appliquee aux points Pi de la 
surface de 1' element rayonnant ER ; 

25 - v 1 correspond au vecteur colonne comportant les valeurs 
des sources S'i de l 1 element rayonnant ER ; et 
v correspond au vecteur colonne comportant les valeurs 
des sources Si de 1* obstacle OBS. 

30 Sur 1' obstacle, la contribution de l'onde emise par 
l f element rayonnant ER s'exprime par : 



31 



V'(P) = P(P).v' t 19 J 

La contribution de l'onde secondaire renvoyee par 

1' Obstacle OBS s'exprime, par definition, par la 
relation : 

V(P) = F(P).v t 2 °l 

Or, dans 1' example represents sur la figure 4A, l'onde 
secondaire correspond simplement a une reflexion de l'onde 
principale. Ce qui s'exprime par la relation : .'• 

V(P) = RV'(P) [21] • 

oil R correspond a une matrice de reflexion dont chaquer 
coefficient represente la contribution a 1' emission, parf 
reflexion, de l'onde secondaire, par chaque source S x (ou : 
Sd, dans le cadre d'une estimation d'une grandeur 
vectorielle) de 1' obstacle OBS. 

Des trois relations [19], [20] et [21], on deduit 
1' expression du vecteur colonne v comportant les valeurs 
des sources de 1' obstacle, a partir du vecteur colonne v 
comportant les valeurs des sources de 1' element rayonnant, 
par la relation : 



v = [F(P)Y l .R.[F'(P)].v' 



[22] 
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Par ailleurs, pour une estimation fine des grandeurs 
scalaires ou vectorielles au point M, notamment pour tenir 
compte de reflexions multiples f il est avantageux de tenir 
compte de la contribution du rayonnement par 1' obstacle, 
au niveau de la surface de l f element rayonnant ER. A cet 
effet, on tient compte, dans 1' estimation des conditions 
aux limites a la surface de 1 T element rayonnant ER (aux 
points P'i) de la contribution du rayonnement des sources 
S'i de 1 T element rayonnant et de la contribution de 
1" emission de l'onde secondaire par les sources Si de 
l 1 obstacle, par la relation : 



On peut ainsi ajuster, gr§ce a la relation [23] , les 
valeurs de sources S'i de 1' element rayonnant ER, en 
tenant compte de la reflexion de 1' obstacle OBS, selon la 
relation suivante : 



Ainsi, on impose simplement des conditions aux limites 
pour l'616ment rayonnant, pour en deduire les valeurs des 
sources v'i* En pratique, on procedera pr6f 6rentiellement 
comme suit : 

- apres maillage des surfaces, on determine la position 
des points Pi et P'i et des sources Si et S'i; 

- en fonction du type d'onde en jeu, on determine les 
coefficients des matrices F(P), F 1 (P) , F(P') et F f (P'); 



[23] 




[10] 
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- en fonction d'une.loi de reflexion de 1' obstacle, on 
determine les coefficients de la matrice de reflexion 
comme dans l'exemple donne plus loin pour une onde 
ultrasonore; 

5 - en fonction de conditions aux liraites sur l 1 element 
rayonnant (dont le comportement est g^neralement connu 
pour un probleme donne) , on determine les valeurs du 
vecteur Vt{P*) aux points P'i de la surface de 1' Element 
rayonnant et on en deduit les valeurs des sources S f ± de. 

10 1' element rayonnant par inversion de la relation [10]; 

- on en deduit aussi les valeurs des sources S± de 
l'obstacle par application de la relation [22]; ■. 

- une fois les valeurs de toutes les sources Si et S'i r ? 
determinees, on peut appliquer le syst&ne matriciel donne > 

15 par la relation [9] a tout point M de l'espace, en.| , 
appliquant a ce point M les matrices d' interaction F et.F 1 ;^ 
(impliquant la position du point. M et les positions des ^ 
sources respectives S± et S'i). L 

20 En se referant a nouveau & la figure 4A, on considere que 
1' obstacle OBS represente simplement une interface entre 
deux milieux Ml et M2, formant ainsi un dioptre qui peut 
§tre plan, tel que represente dans l'exemple de la figure 
4A, mais aussi incurve ou de forme generale quelconque. 

25 Les coefficients de reflexion R± associes a chaque point 
Pi dependent, dans le cadre de la propagation d'une onde 
ultrasonore ou electromagnetique de haute frequence, de 
I 1 angle d T incidence Pi du rayon issu de la source Si, au 
point de l f espace tridimensionnel M. 

30 
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Pour une onde ultrasonore, 1' expression des coefficients 
de reflexion Ri est donnee par : 



R; = 



p 2 c 2 cosy?/ - p\C\ 



1 _4 + 4cos 2 /3 / 



1/2 



p 2 C 2 COSyff/ + P\C\ 



1 _4 + 4cos 2 A - 



1/2 



[24] 



ou : 

- ci est la vitesse du son dans le milieu Mi ; 

- C2 est la vitesse du son dans le milieu M 2 ; 

- pi est la densite du milieu M a ; 
10 - p2 est la densite du milieu M 2 . 

Dans cette expression [24], le terme cospi peut simplement 
etre estime en fonction des coordonnees dans l'espace du 
point M et du point representant la source Si. 

15 En se referant maintenant a la figure 4B, la meme 
estimation peut etre menee pour un point M situe dans le 
milieu M 2 . Dans ce cas, l'onde que recoit le point M est 
une onde transmise par 1' obstacle OBS. En particulier, on 
remarque que les sources de 1' element rayonnant ER ne sont 

20 plus actives, du fait de 1 ' occultation de 1' element 
rayonnant ER par 1' obstacle OBS. En transmission, le 
raisonnement s' applique coitime precedemment avec une 
condition aux limites imposee aux points Pi par les 
valeurs des coefficients de transmission Tj. associes a 

25 chaque point Pi- Dans le cadre de la propagation d'une 
onde ultrasonore, chaque coefficient de transmission Ti 
est donne par la relation : 
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2/? 2 c 2 cosfl [25] 

lj — 







2 


1/2 


p 2 c 2 COSjffj + pici 




-§• cos 2 Pi 













Comme indique precedemment, les termes cospi peuvent etre 
determines en fonction des coordonnees respectives des 
sources Si et du point M. 

Pour estimer les valeurs de sources Si de 1' obstacle OBS, 
on applique la " relation [22] en remplacant toutefois la 
matrice de reflexion R par la matrice de transmission T : -'r 

» 

v^fFCP^TtPCP)]^ [26] t 

Dans le cadre de la propagation d'une onde ultrasonore,j 
les coefficients des matrices R et T s'estiment pour'' 
chaque source S ± et pour chaque point Pi- En particulier, 
chaque coefficient T ifj ou R ifj de la matrice T ou de la 
matrice R (ou i correspond a la i feme ligne et j correspond 
a la j 6me colonne) s'exprime en fonction d'un angle fa 
entre une normale a la surface de 1' obstacle au point P± 
et une droite passant par le point Pi et par une source 
S 3 . On peut ecrire ainsi, de maniere generale, les deux 
relations exprimant les valeurs des coefficients des 
matrices T et R par les relations respectives suivantes : 

Ty = f t (cos/?ij) [27] 
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Rij = f r (cos/?ij) 



[28] 



ou f t est donne par la relation [25] et f r est donne par 
la relation [24] . 

5 

De fa<?on plus generale, en reference aux figures 4A et 4B, 
on indique que, si l'on considdre 1" obstacle comme un 
materiau plein representant un milieu M2 distinct d'un 
milieu Ml dans lequel se propageait initialement l'onde 

10 principale : 

- pour une reflexion de l'onde principale sur 1 T obstacle 
en tant que milieu M2 (la surface de l 1 obstacle formant 
dioptre entre les milieux Ml et M2) f les hemispheres HEMi 
sont orientes vers l'exterieur de l 1 obstacle (figure 4A) ; 

15 - pour une transmission de l'onde principale dans 
l 1 obstacle, les hemispheres HEMi sont orientes vers 
l'interieur de 1' obstacle (figure 4B) . 

On se refere maintenant a la figure 5A pour decrire le cas 
20 d'un 1' obstacle plan OBS de dimensions finies, excite par 
un element rayonnant ER, incline d f un angle predetermine 
par rapport a 1' obstacle OBS. Comme indiqu<§ precedemment , 
pour une onde ultrasonore, on prendra en compte 
1 1 inclinaison de 1' Element rayonnant pour calculer la 
25 contribution de l'onde 6mise par l f element rayonnant au 
point M. Par ailleurs, de fagon particulierement 
avantageuse, on maille une surface qui englobe la surface 
de 1 T obstacle (figure 5A) . Pour une tranche de l'espace 
delimit^e par l f element rayonnant, d'une part, et 
30 l'obstacle, d'autre part (figure 5A) , on peut consid^rer 
trois types de sources : 
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les sources S'i de 1' element rayonnant ER, 
des sources S0 if qui renvoient l'onde secondaire, par 
reflexion de 1' obstacle OBS, en fonction d'un certain 
coefficient de reflexion R de 1' obstacle ; et 
des sources SSi, qui ne renvoient pas d'onde secondaire 
et auxquelles un coefficient de reflexion nul peut etre 
attribue si 1' obstacle separe deux milieux d' indices 
identiques. Dans ce cas, ces sources SSi sont 
considerees comme "eteintes" dans la tranche de 
l'espace precitee et ne sont pas prises en compte dans 
les calculs de la grandeur physique au point M de la 
figure 5A. En revanche, ces sources SSi peuvent §tre^ 
actives par reflexion de l'onde principale si^ 
1' obstacle OBS separe deux milieux d' indices^ 
differents. 

Par ailleurs, pour estimer les grandeurs scalaires otf 
vectorielles associees a un point M d'un demi-espace 
delimite par la surface englobant 1' obstacle OBS (a droite 
de la figure 5B) , on considere : 

- les sources SO'i.de l'obstacle, actives par transmission 
de l'onde principale, et 

- les sources SSi, auxquelles on affecte maintenant un 
coefficient de transmission egal a 1 si 1' obstacle 
separe deux milieux de memes indices. Ces sources SSi se 
comportent finalement (aux angles d' incidence pres) 
comme les sources S'i de 1 ' element rayonnant ER. 

Les sources S'i de 1' element rayonnant peuvent alors etre 
"eteintes" pour le calcul des grandeurs physiques dans ce 
demi-espace. 
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Pour calculer les valeurs v r des sources S'i de l x Element 
rayonnant (desquelles sont deduites les valeurs v des 
sources de l'obstacle selon les relations [22] et [26]), 
on appliquera s implement des conditions aux limites aux 
5 points de la surface active de l 1 element rayonnant ER. Par 
exemple f pour une propagation d'onde ultrasonore, on peut 
indiquer que les vitesses sonores aux points de la surface 
de 1' Element rayonnant ER sont perpendiculaires a cette 
surface et de modules v 0 egaux entre eux. 

10 

De mani^re g6nerale, on indique que l'espace 
tridimensionnel peut etre ainsi decoupe par des interfaces 
d61imitant des milieux de proprietes distinctes, chaque 
interface representant un obstacle au sens de la pr§sente 

15 invention. On calcule alors les grandeurs physiques dans 
chaque tranche de I'espace. Par exemple, dans le cadre de 
1' etude d f uhe h£terostructure (a plusieurs interfaces) on 
applique le procede ci-avant pour des tranches successives 
de l'espace en considerant deux interfaces : l'une 

20 representant un "element rayonnant" au sens des figures 4A 
et 5A, par exemple par transmission d'une onde regue, et 
1' autre representant un obstacle recevant I'onde 
transmise. Avantageusement, on tient compte, pour chaque 
tranche de l'espace, des contributions de toutes les 

25 interfaces comme exprime par les relations [10] et [22] . 

La simulation de la figure 5C correspond, pour une onde 
ultrasonore, k la situation des figures 5A et 5B en tenant 
compte : 

30 - de la contribution de l f emission de l f onde principale 
par l 1 element rayonnant ER ; 
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- de la contribution de la reflexion de cette onde 
principale par 1' obstacle ; et 

- de la contribution de la transmission de l'onde 
principale par 1' obstacle. 

Les lignes de niveau de la figure 5C correspondent a 
differents paliers de pression acoustique. L» element 
rayonnant ER est dispose a 10 mm. de 1' obstacle. OBS et 
incline de 20° par rapport a ce dernier. On remarque en 
particulier des f ranges d' interferences dans une zone 
entre 1' obstacle OBS et 1' element rayonnant ER. Une telle 
simulation peut avantageusement indiquer une position, 
ideale d'un capteur ultrasonore. Ces capteurs ultrasonores 
.comportent habituellement un transducteur comme element 
rayonnant actif et un detecteur pour mesurer les onde^ 
ultrasonores recues. La simulation de la figure 5C peu^ 
ainsi indiquer en outre la forme ideale d'un capteur. 
ultrasonore, selon les applications souhaitees, pour une. 
forme d' obstacle donnee. 



La simulation de la figure 5C a ete effectuee grace a un 
calcul matriciel programme a l'aide du logiciel de calcul 
MATLAB© . Le nombre .. de mailles total choisi pour 
l'obstacle et pour 1' element rayonnant (ici, quelques 
25 centaines en tout) est alors optimise : • 

- d'une part, pour limiter la duree des calculs ; et 

- d' autre part, pour que la taille des mailles reste 
inferieure a une demie longueur d'onde, de maniere a 
verifier le critere de Rayleigh. 
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On indique toutefois que, comme les elements a mailler 
dans la mise en ceuvre du procede selon 1 1 invention sont 
s implement des surfaces, les temps de calcul sont bien 
moins longs que ceux necessaires dans la mise en oeuvre 
5 d'un procede de calcul de type "par elements finis". 

La presente invention peut ainsi se manifester par la mise 
en oeuvre d'une succession d 1 instructions d'un produit 
programme d'ordinateur stocks dans la memoire d'un disque 
10 dur ou d'un support amovible et se deroulant comme suit : 

- choix d'un pas de maillage notamment en fonction de la 
longueur d'onde de l'onde principale ; 

- determination des coordonnees des sources S± et/ou S'i 
et des points de contact Pi et/ou P'i ; 

15 - choix d ' un type d ' onde en j eu et calcul des 
coefficients des matrices d' interaction appliquees aux 
points Pi et/ou P'i par la mise en oeuvre d'un logiciel 
de calcul matriciel ; 

- choix d'une loi de reflexion et/ou de transmission de 
20 la surface de l f obstacle et calcul des coefficients des 

matrices de reflexion et/ou de transmission ; 

- calcul des valeurs des sources Si et/ou S'± ; et 

- calcul des grandeurs scalaires ou vectorielles en tout 
point de 1 ' espace t ridimensionnel <. 

25 

A ce titre, la presente invention vise aussi un tel 
produit programme d'ordinateur. 

Bien entendu, la presente invention ne se limite pas £ la 
30 forme de realisation d^crite ci-avant a titre d'exemple ; 
elle s'etend & d'autres variantes- 
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Ainsi, on comprendra que, meme si dans les figures 
commentees ci-avant . on represente a la fois la surface 
d'un obstacle et la surface d'un element rayonnant, la 
presente invention s' applique aussi a 1 ' estimation de 
grandeurs physiques dans le cadre d'une onde interagissant 
avec un obstacle et emise en champ lointain. Dans ce 
contexte, il n'est pas necessaire de materialiser la 
surface d'un element rayonnant a mailler et les relations 
[8] et [14] ci-avant suffisent a determiner 1' interaction 
entre cette onde et 1' obstacle. 
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On a indique ci-avant des equations permettant de calculer 
les grandeurs - scalaires ou vectorielles en un point M dj 
l'espace, pour des ondes electromagnetiques, ou encore 
acoustiques. Bien- entendu, ces grandeurs peuvent §tre 
estimees pour d'autres types d' ondes, notamment pour des 
ondes thermiques, des ondes electromagnetiques impliquani: 
des antennes radiof requences, ou autres. 

Bien .entendu, la presente invention ne se limite pas a une 
application au contrdle non destructif, mais a tout type 
d' application, notamment en imagerie medicale, par exemple 
pour 1 ' 6tude de microsystemes mettant en oeuvre une 
microscopie acoustique a miroirs mobiles. 

On a decrit ci-avant des interactions entre une onde et un 
obstacle unique. 'Bien entendu, la presente invention 
s' applique a une interaction avec plusieurs obstacles. A 
cet effet, il convient simplement de mailler les surfaces 
de ces obstacles et d'additionner leur contribution pour 
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1* estimation d'une grandeur vectorielle ou scalaire en un 
point quelconque de l'espace. De meme, comme indique ci- 
avant, la surface de 1' obstacle OBS peut etre plane, ou 
encore incurv^e, ou encore de forme complexe quelconque. 

5 

Ainsi, dans le cadre d'une onde interagissant avec 
plusieurs obstacles dans l'espace, une simulation 
equivalente a celle representee sur la figure 5C 
permettrait de positionner des capteurs et/ou des elements 
10 rayonnants en fonction de la configuration de ces 
obstacles, notamment pour une application a la 
determination de la position de hauts-parleurs dans un 
habitacle cloisonn§, tel qu'un habitacle de vehicule 
automobile ♦ 

15 

L'espace tridimensionnel peut etre decoupe en une 
pluralite de regions, comme decrit ci-avant en reference 
aux figures 4A, 4B, 5A et 5B. Toutefois, pour qu'une 
surface de 1'une desdites regions soit consideree comme un 
20 obstacle d'une onde principale, devenant actif par 
Emission d'une onde secondaire, 1' incidence de 1 1 onde 
principale sur cette surface doit pref erentiellement 
rester inf^rieure ou 6gale a 90°. 
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Revendi cations 

1* Procede pour evaluer une grandeur physique associee a 
une interaction entre une onde et un obstacle , dans une 
region de I'espace tridimensionnel, dans lequel : 

a) on determine, par maillage, une plurality 
d'^chantill.ons de surface (dSi) dont une partie. au moins 
represente la surface d'un obstacle recevant une onde 
principale et emettant, en reponse, une onde secondaire, 
et I'on attribue a chaque echantillon de surface au moins 
une source (Si) emettant une onde elementaire representant 
une contribution a ladite onde secondaire, A 

b) on forme un systeme matriciel comportant : 

- une matrice d' interaction (F(M)), inversible^ 
appliquee a une region donnee (M) de I'espace et 
comportant un nombre de colonnes correspondant a ua 
nombre total de sources, .£ 
une premiere matrice colonne dont chaque coefficient 
(vi) est associe* a une source (Si) et caracterise 
l f onde elementaire qu T elle <§met, 

- et une seconde matrice colonne, obtenue par une 
multiplication de la premiere matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d'une grandeur physique (V(M)) 
representative de l'onde emise par l f ensemble des 
sources en ladite region donnee (M) , 

c) pour estimer les coefficients de la premiere matrice 
colonne (Vi) , on affecte des valeurs de grandeur physique 
choisies (V(Pi)) a des points predetermines (Pi), associes 
chacun d un echantillon de surface (dS ± ) , lesdites valeurs 
choisies (V(Pi)) 6tant rangees dans la seconde matrice 
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colonne, et l'on multiplie cette seconde matrice colonne 
par l 1 inverse de la matrice d 1 interaction appliquee 
auxdits points predetermines (Pi), 

d) pour evaluer ladite grandeur physique (V(M)) 
5 representant l'onde emise par 1' ensemble des sources dans 
une region donnee (M) de l'espace tridimensional, on 
applique la matrice d' interaction a ladite region donnee 
<M) et on multiplie cette matrice d' interaction par la 
premiere matrice colonne comportant les coefficients 
10 estim6s a l'etape c) . 

2. Procede selon la revendication 1, dans lequel, pour 
6valuer une grandeur physique representative d'une 
interaction entre un element rayonnant une onde principale 
15 et un obstacle recevant cette onde principale, 

- & l'6tape a) , on determine en outre, par maillage, une 
pluralite d' echantillons de surface (dS f i) representant 
ensemble une surface active de 1' element rayonnant l'onde 
principale, et l'on attribue a chaque echantillon de la 

20 surface active au moins une source (S'i) emettant une onde 
elementaire representant une contribution a ladite onde 
principale, 

- on applique en outre les etapes b) , c) et d) aux 
echantillons de la surface active, et 

25 - on lvalue ladite grandeur physique (V(M) ) representant 
l 1 interaction entre l'616ment rayonnant et I 1 obstacle dans 
une region donn6e (M) de l'espace tridimensionnel , en 
prenant en compte la contribution, en ladite region donn§e 
(M) , de l'onde principale emise par l 1 ensemble des sources 

30 de la surface active et la contribution de l'onde 
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secondaire emise par 1' ensemble des sources de la surface 
de 1' obstacle. 

3. Procede selon l'une des revendications 1 et 2, dans 
lequel chaque coefficient de la matrice d' interaction, 
appliquee a une region donnee de l'espace, est 
representatif d'une interaction entre une source et ladite 
region donnee et la . valeur de chaque coefficient est 
fonction d'une distance entre une source et ladite region 
donnee . 

4. Procede selon l'une des revendications 1 a 3, dans 
lequel la matrice d' interaction appliquee", a l'etape ok, 
auxdits points predetermines (Pi), comporte un nombre de 
lignes correspondant a un. nombre total de points 
predetermines (Pi) . & 

5. Procede selon l'une des revendications 1 a 4, dar^s 
lequel la grandeur physique a evaluer est une grandeur 
scalaire (V(Pi)) et, a l'.etape a), on attribue une source 
unique a chaque echantillon de surface. 

6. Procede selon la revendication 5, dans lequel la 
matrice d' interaction (F (M) ) appliquee, a l'etape d) , a 
une region de l'espace (M) , comporte une ligne. 

7. Procede selon l'une des revendications 5 et 6, dans 
lequel chaque point predetermine (Pi) associe a un 
.echantillon de surface (dSi) correspond a un point de 

• contact entre cet echantillon de surface (dSi) et un 
hemisphere dont la surface est egale a la surface de cet 
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echantillon de surface, et de centre correspondant a une 
position de la source (Si) qui est attribute & cet 
echantillon de surface. 

5 8. Precede selon l'une des revendications 5 a 7, dans 
lequel : 

- l'onde principale est une onde electrique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge electrique associees chacune k une 

10 source, et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel electrique. 

9. Procede selon l'une des revendications 5 a 7, dans 
15 lequel : 

- l'onde principale est une onde magnetique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de flux magnetique associees chacune a une source, 
et 

20 - les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel magnetique. 

10. Procede selon l'une des revendications 5a 7, dans 
lequel: 

25 - l'onde principale est une onde sonore, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de vitesse de son associees chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
30 valeurs de pression acoustique. 
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11. Proceed selon l'une des revendications 1 a 4, dans 
lequel la grandeur physique a evaluer est une grandeur 
vectorielle (V(P±)) exprimee par ses trois coordonnees 
dans l'espace tridimensionnel, et l'on attribue, a l'etape 
a), trois sources (SAi, SBi, SC A ) a chaque echantillon de 
surface (dSi) . 

12. Procede selon la revendication 11, dans lequel la 
matrice d' interaction (F V (M)) appliquee, a l'etape d) , k 
une region de l'espace (M) , comporte une ligne pour chaque 
coordonnee <X,Y,Z) de l'espace. 

13. Procede selon l'une des revendications 11 et 12, dans 

lequel : h 
-les trois sources attributes a chaque echantillon.de 

surface sont sensiblement dans un meme plan et ^ 
- chaque point predetermine (Pi) associe a un echantillon 
de surface (dSi) correspond a un point de contact entrj 
cet echantillon et un hemisphere dont la surface est egale 
a la surface de . cet echantillon, et de centre 
correspondant a la position d'un barycentre des trois 
sources . 

14. Procede selon la revendication 13, dans lequel les 
trois sources d'un meme Echantillon de surface forment 
sensiblement un triangle Equilateral, et les triangles des 
Ecnantillons de surface sont orientes sensiblement 
aleatoirement les uns par rapport aux autres. 

15. Procede selon l'une des revendications 11 a 14, dans 
lequel : 
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- l'onde principale est une onde elect rique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge electrique associees chacune a une 
source, et 

5 - les coefficients de la seconde matrice colonne sont des 
valeurs de coordonnees d f un champ electrique. 

16. Procede selon l'une des revendications 11 a 14 , dans 
lequel : 

10 - I'onde principale est une onde magnetique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de flux magnetique associees chacune a une source,, 
et 

- les coefficients de la seconde matrice colonne sont des 
15 valeurs de coordonnees d'un champ magnetique. 

17. Proced6 selon l'une des revendications 11 a 14 , dans 
lequel : 

- l'onde principale est une onde sonore, 

20 - les coefficients de la premiere matrice colonne sont des 
valeurs de vitesse de son associees chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de coordonnees d'une vitesse acoustique. 

25 

18. Procede selon l'une des revendications pr£cedentes, 
dans lequel, pour estimer la contribution de l'onde 
secondaire en ladite region donnee a l'etape d) , lesdites 
valeurs de grandeur physique (V(P A )) choisies a l'etape c) 

30 sont fonction d'un coefficient predetermine de reflexion 
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et/ou de transmission de l'onde principale par chaque 
echantillon de surface de 1' obstacle. 

19. Procede selon la revendication 18, prise en 
combinaison avec l'une des revendications 6 et 12, dans 
lequel l'onde secondaire correspond a une reflexion de 
l'onde principale sur 1' obstacle et 1 ' hemisphere est 
oriente vers l'exterieur de 1.' obstacle. 

20. Procede selon la revendication 18, prise en 
combinaison avec l'une des revendications 6 et 12, dans 
lequel l'onde secondaire correspond a une transmission de 
l'onde principale dans 1' obstacle et 1' hemisphere est 
oriente vers l'interieur de 1' obstacle. 

21. Procede selon l'une des revendications 19 et 20, dans 
lequel, a l'etape c) , on determine les valeurs (v'i) 
associees aux sources (S' ± ) de 1' element rayonnant (ER) e£ 
on met en forme au moins : 

- une premiere matrice d' interaction (F(P)) representant 
la contribution des sources de 1' obstacle aux points 
predetermines de la surface de 1' obstacle (Pi), 

- une seconde matrice d' interaction (F' (P) ) representant 
la contribution . des sources de 1' element rayonnant aux 
points predetermines de la surface de 1' obstacle (Pi), 

- une matrice de reflexion (R) ou de transmission (T) , 
dont les coefficients representent des coefficients de 
reflexion ou de transmission en chaque point predetermine 
(P ± ) de 1' obstacle, 

pour determiner les valeurs des sources de 1' obstacle (vi) 
en fonction des valeurs des sources de 1» element rayonnant 
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(v'i) et d'une multiplication des premiere et seconde 
matrices d' interaction et de la matrice de reflexion ou de 
transmission . 

5 22. Procede selon la revindication 21, dans lequel, a 
l'etape c) , on determine les valeurs (v'i) associees aux 
sources (S'i) de l 1 element rayonnant (ER) en tenant compte 
de la reception de I'onde secondaire par 1' Element 
rayonnant (ER) et en mettant en forme en outre : 

10 - une troisieme matrice d T interaction (F(P')) representant 
la contribution des sources de l 1 obstacle aux points 
predetermines de la surface de 1' element rayonnant (P'i), 

et une quatrieme matrice d 1 interaction (F' (P T ) ) 
representant la contribution des sources de 1 T element 

15 rayonnant aux points predetermines de la surface de 
1* element rayonnant (P'i) . 

23. Precede selon 1'une des revendications 19 a 22, dans 
lequel la surface de l 1 obstacle correspond a une interface 

20 entre deux milieux distincts d'une heterostructure. 

24. Procede selon 1'une des revendications precedentes, 
dans lequel I'onde principale est une onde sonore et les 
coefficients de la matrice d' interaction sont chacun 

25 fonction d'un angle d' incidence d'une onde eiementaire 
issue d'une source en ladite region donnee (M) . 

25. Procede selon 1'une des revendications precedentes, 
dans lequel l f onde principale est une onde sonore et, a 

30 l'etape a), on choisit un nombre total d'echantilions de 
surface (dSi) sensiblement en fonction d'une longueur 
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d'onde de l'onde sonore pour verifier le critere de 
Rayleigh. 

26. Procede selon l'une des revendications precedentes, 
dans lequel on compare une pluralite de valeurs de la 
grandeur physique estimee a l'etape d) , obtenues pour une 
pluralite de regions de l'espace, pour selectionner une 
region candidate pour la disposition d'un element 
• rayonnant destine a interagir avec 1' obstacle. 



27. Procede selon l'une des revendications 2 a 26, dans 

lequel 1' element rayonnant est un capteur, pour un 

contrSle non destructif, destine a analyser un obje|: 

formant un . obstacle de l'onde principale. . ± 

15 ir- 
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